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Abstract 

We study the stochastic dynamics of c and b quarks in the hot plasma produced in nucleus-nucleus collisions at RHIC 
and LHC, providing results for the nuclear modification factor Raa and the elliptic flow coefficient V2 of the single- 
electron spectra arising from their semi-leptonic decays. The initial QQ pairs are generated using the POWHEG 
code, implementing pQCD at NLO. For the propagation in the plasma we develop a relativistic Langevin equation 
(solved in a medium described by hydrodynamics) whose transport coefficients are evaluated through a first-principle 
calculation. Finally, at T c , the heavy quarks are made hadronize and decay into electrons: the resulting spectra are 
then compared with RHIC results. Predictions for LHC are also attempted. 

Keywords: Quark Gluon Plasma, heavy quarks, Langevin equation, transport coefficients 



1. Introduction 

Heavy quarks, produced in initial hard processes, allow to perform a "tomography" of the medium created in the 
relativistic heavy-ion collisions at RHIC and (soon) at LHC. The modification of their spectra provides information on 
the properties (encoded into few transport coefficients) of the matter (hopefully a thermalized Quark Gluon Plasma) 
crossed before hadronizing and giving rise to experimental signals: so far, the electrons from their semi-leptonic 
decays. We employ an approach based on the relativistic Langevin equation, assuming that medium-modifications of 
the initial heavy-quark spectrum arise from the effects of many independent random collisions. As a final outcome we 
provide results for the Raa and v'2 of non-photonic electrons from heavy-flavor decays measured at RHIC and attempt 
predictions for LHC. Further results can be found in Refs. JH01- For similar studies see Refs. JH 0, El S 01 • 

2. The relativistic Langevin equation 

The usual Langevin equation can be generalized to the relativistic case 112], providing a tool to study the propa- 
gation of c and b quarks in the QGP. The variation of the heavy-quark momentum in the time-interval At 

^t=-r lD {p)p i +af), (1) 

is given by the sum of a deterministic friction term and a stochastic noise term %'(t), which is completely determined 
by its two-point temporal correlator 



(?(t)¥(t)) = b'j(p)8(t - t'\ with b ij (p) = K L (p)p i p j + K T (p)(5 ij - p l p j ). (2) 

The latter involves the transport coefficients Kj{p) = and Ki(p) = ^j^, representing the average transverse 

and longitudinal squared-momentum acquired per unit time by the heavy quark due to the collisions suffered in the 
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Figure 1: The transport coefficients Kt/l of c and b quarks after summing the soft and hard contributions. The dependence on the cutoff |?|* is mild. 
The coupling g was evaluated at the "central value" of our systematic scan /i = 1 .5nT . 



medium. Finally, as in the non-relativistic case, the friction coefficient rjoip) is fixed in order to insure the approach 
to thermal equilibrium. In the Ito discretization Jit] of Eq. ([TJ one has: 
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Eq. (Q]) has then to be solved in the evolving medium produced in the heavy-ion collisions and described by ideal/viscous 
hydrodynamics. For this purpose the output of two independent hydro codes SEED] are exploited. Details on the 
employed procedure can be found in Refs. fl 



3. Evaluation of the transport coefficients 

The transport coefficients kt/l are evaluated according to the procedure presented in Refs. (jl 0]. Following 
Ref. lfl3ll we introduce an intermediate cutoff \t\* ~ m 2 D (t = (P' - P) 2 ) separating hard and soft scatterings. The 
contribution of hard collisions (|f| > \t\*) is evaluated through a kinetic pQCD calculation of the processes Q(P)q i /- i — > 
Q(P')qi/i and Q(P)g -> Q(P')g u3- On the other hand in soft collisions (|f| < \t\*) the exchanged gluon feels the 
presence of the plasma. A resummation of medium effects is thus required and this is provided by the Hard Thermal 
Loop approximation. The final result is given by the sum of the two contributions kt/l(p) = K ^pt(p) + *ryx(p) an d i ts 
explicit expression can be found in Ref. In Fig.Q]we display the behavior of the transport coefficients of c and b 
quarks. The sensitivity to the value of the intermediate cutoff |f|* is quite small, hence supporting the validity of the 
approach. 



4. Numerical results: single-electron spectra at RH1C and LHC 

For each explored case we generated an initial sample of 45 • 10 6 cc and bb pairs, using the POWHEG code 111 511 
with CTEQ6M PDFs. In the AA case we introduced nuclear effects in the PDFs according to the EPS09 scheme [16]; 
the quarks were then distributed in the transverse plane according to the nuclear overlap function dN/dx± ~ Tab(x, y) = 
TA(x + b/2,y)TB(x-b/2,y). At the proper-time r = t 2 -z 2 = tq we started following the Langevin dynamics of the 
quarks until hadronization. The latter was modeled using Peterson fragmentation functions lfl7ll . with branching 
fractions into the different hadrons taken from Refs. lfl8i[l9ll . Finally each hadron was forced to decay into electrons 



with PYTHIA [20], using updated decay tables 12111 , The e-spectra from c and b were then combined with a weight 



accounting for the respective total production cross-section and for the branching ratios of the various processes 
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Table 1: Initialization for RHIC and LHC: the cc and bb production cross section given by POWHEG and the explored hydro scenarios. 
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Figure 2: The electron R^i for minimum bias Au-Au collisions (corresponding to an impact parameter i> = 8.44 fm) at V*NN =200 GeV compared 
with PHENIX results 1231 . The sensitivity to the scale n=nT - 2nT is displayed. The high-momentum region (pj >4 GeV/c) is better reproduced 
with the intermediate coupling. Ideal (left panel) and viscous (right panel) hydro scenarios provide similar results. 




Figure 3: The electron Raa and vi for minimum bias Pb-Pb collisions (b = 8.58 fm) at LHC for the two hydro scenarios considered. 



leading to electrons in the final state. We first address the RHIC case. The parameters characterizing the "initial 
state" are given in Table Q] In Fig. [2] our results for the electron Raa in minimum-bias collisions are displayed and 
compared with PHENIX data B22I1 . The dependence on the hydro scenario is quite small. The major uncertainty 
comes from the scale at which one evaluates the coupling g: for that we explore the range ji-nT - 2nT . The high- 
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momentum region (pj > 4 GeV/c) is better reproduced evaluating a s at the scale 1 .5nT (a s — 0.32 at T — 300 Me V). 
Hadronization via coalescence [4], so far ignored, could improve the agreement at lower pj- Results for v% can be 
found in Ref. yfl. We address now the LHC case, described by the parameters in Table Q] Notice the dramatic effect 
of the nPDFs on the cc production. Our findings for Raa and \>2 are shown in Fig. [3] While the Z?-spectrum is simply 
suppressed at high-p^, the charm - beside a stronger quenching - also displays a sizable flow. Finally in Fig. |4]we 
provide a comparison of our results for the inclusive (c + b) single-electron spectra for RHIC and LHC conditions. In 
the last case the larger initial temperature entails a stronger quenching and a more pronounced elliptic flow. The effect 
would be even stronger in the absence of initial state effects, responsible for a larger cr^/cr cf at LHC with respect to 
RHIC. 
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